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transformed path measure under a trap Hamiltonian, which prevents the motions from 
escaping to infinity, and a pair-interaction Hamiltonian, which imposes a repellency of 
the N paths. In fact, this interaction is an A-dependent regularisation of the Brown¬ 
ian intersection local times, an object which is of independent interest in the theory of 
stochastic processes. 

The time horizon (interpreted as the inverse temperature) is kept fixed. We analyse 
the model for diverging number of Brownian motions in terms of a large deviation 
principle. The resulting variational formula is the positive-temperature analogue of the 
well-known Gross-Pitaevskii formula, which approximates the ground state of a certain 
dilute large quantum system; the kinetic energy term of that formula is replaced by a 
probabilistic energy functional. 

This study is a continuation of the analysis in [ABKOdj where we considered the 
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MSC 2000. 60F10; 60J65; 82B10; 82B26. 

Keywords and phrases. Interacting Brownian motions, Brownian intersection local times, large devi¬ 
ations, occupation measure, Gross-Pitaevskii formula. 

^Max-Planck Institute for Mathematics in the Sciences, Inselstrafie 22-26, D-04103 Leipzig, Germany, 
adamsSmis.mpg.de 

^Fachbereich Mathematik und Informatik, Johannes-Gutenberg-Universitat Mainz, Staudingerweg 9, D-55099 Mainz, 
Germany, jbbruSmathematik. uni-mainz. de 

^Mathematisches Institut, Universitat Leipzig, Augustusplatz 10/11, D-04109 Leipzig, Germany, 

koenigOmath.uni-leipzig.de 





2 


STEFAN ADAMS, JEAN-BERNARD BRU AND WOLFGANG KONIG 


1. Introduction and main results 


1.1 Background. 

One of the important problems in mathematical physics is the understanding of large systems of 
interacting quantum particles at extremely low or zero temperature. This question is raised in the 
literature since decades and is studied under a lot of different view points. Large systems of bosons (i.e., 
quantum particle systems whose wave functions are invariant under permutation of the single-particle 
variables), at extremely low temperatures, may undergo a phase transition, the so-called Bose-Einstein 
eondensation: a macroscopic portion of the system is described by one suitable single-particle wave 
function. This phenomenon was shown for ideal (i.e., non-interacting) gases already in 1925 by Bose 
and Einstein. A rigorous understanding of condensation effects in various ultra cold materials remains 
a widely open and challenging problem until today. This question received an enormous impetus by 
the first experimental realisation of Bose-Einstein condensation in 1995. 

A many-particle quantum system is mathematically described by a A^-particle Hamilton operator 
consisting of a kinetic energy term, a trap term and an interaction term. Its spectral analysis, at 
least for realistic interacting models, is out of reach of contemporary analysis. Rigorous theoretical 
research started with Bogoliubov and Landau in the 1940ies, followed by Penrose, Eeynman and many 
others. They analysed simplified mathematical models featuring only the most important physical 
phenomena. However, these approaches turned out to be intuitively appealing and relevant. See 
jABOdalfTBlMhl for a review and some recent results. 

Another mathematical approach is to consider systems that are dilute on a particular scale and 
are kept within a bounded region by the presence of a trap. Here ‘dilute’ means that the range 
of the interparticle interaction is small compared with the mean particle distance. These systems 
are supposed to be easier to analyse at least as it concerns the ground state. In a particular dilute 
situation, the ground states and their energy were analysed in the many-particle limit |IjSSYn5] . It 
turned out that the well-known Gross-Pitaevskii formula describes the system remarkably well. This 
variational formula has a kinetic term (the usual energy), a trap term and a quartic term with a 
pre factor. As was predicted by earlier theoretical work, the only parameter of the pair interaction 
functional that persists in the limit is its scattering length. See |PSn3j for an overview about the 
physics and |IjSSYn5] for an account on recent mathematical research. 

However, the mathematically rigorous understanding of large quantum systems at positive temper¬ 
ature is still incomplete. For dilute systems of fermions (i.e., quantum particle systems whose wave 
functions are antisymmetric under permutation of the single-particle variables), first results for posi¬ 
tive temperature are in [Sein5j . One main concern of quantum statistical mechanics is to evaluate the 
trace of the Boltzmann factor for inverse temperature /3 > 0 to calculate all thermodynamic 

functions. The Feynman-Kae formula provides a representation of these traces as functional integrals 
over the snace of Browniau naths on the finite time horizon [0. B] [din m. Hence, it is clear that an 
appropriate description of quantum systems at positive temperature is given in terms of independent 
Brownian motions in a trap with a mutually repellent pair interaction. The trap and the interaction 
are imposed via exponential densities, so-called Hamiltonians. 

In the present paper, we make a contribution to a rigorous analysis of a certain model of a large 
number of mutually repellent Brownian motions in a trap at any positive temperature. We introduced 
this model in earlier work |ABKn4j . The pair interaction in that model is a path interaction, not a 
particle interaction. It turned out there that its behaviour in the zero-temperature limit is asymptot¬ 
ically well described by a variational formula known as the Hartree formula. Therefore, we call this 
Brownian model the Hartree model. The interaction Hamiltonian is given via a double time integra¬ 
tion and thus the Hartree model is related to Polaron type models |DV83) , |BDS98] , where instead of 
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several paths a single path is considered. In [ARKOdj we showed that the many-particle limit of the 
Hartree formula is well approximated by the above mentioned Gross-Pitaevskii formula. However, the 
decisive parameter here is not the scattering length, but the integral of the pair interaction functional. 

While that result describes the zero-temperature situation, in the present work we study the case 
of positive temperature, i.e., interacting Brownian motions on a fixed finite time horizon. Our main 
result is a description of the many-particle limit in terms of a variational formula that is analogous to 
the Gross-Pitaevskii formula. The only difference is the kinetic energy term, which is now replaced by 
a probabilistic energy term. In the language of large deviation theory, this term is the rate function 
that governs a certain large deviation principle. In our particular case it is a certain Legendre-Fenchel 
transform of an exponential Brownian integral. 

The system of motions is dilute on the same scale as the above mentioned models. This means in 
our particular case that the repellent interaction is an approximation of a highly irregular functional 
of the motions, the so-called Brownian intersection local times, which measure the amount of time 
that is spent by two motions at their intersection points. This object is of independent interest in the 
theory of stochastic processes. Our main result implicitly states a large deviation principle for the 
mean of regularisations of the intersection local times, taken over all mutual intersections of a large 
number of Brownian paths. 

The remainder of Section ^ is organised as follows. We introduce the model in Section 11.21 and 
present our main result and some conclusions in Section fl.dl In Section El we embed these results in 
a broader perspective, discuss our results and mention some open problems. 

The remainder of the paper is structured as follows. In Section H we prove some properties of the 
probabilistic energy term and the variational formula. Section H contains the proof of our main result. 
In the Appendix we give a short account on large deviation theory in Section f4.II and recall a related 
result by Lieb et al. on the large-limit of the ground state in Section 14.21 

1.2 The model. 

We consider a family of N independent Brownian motions, ..., i> 0 , in with gener¬ 

ator —A each. We assume that each motion possesses the same initial distribution, which we do not 
want to make explicit. The model we study is defined in terms of a Hamiltonian which consists of two 
parts: a trap part, 

A 

= WiBi'^)ds, ( 1 . 1 ) 

i=i -^0 

and a pair-interaction part, 

= E ^ r r (1-2) 

Jo 

Here IT: —> [0, oo] is the so-called trap potential satisfying {x) = oo, and v. (0,oo) —> 

[0,oo] is a pair-interaction function satisfying 0 < limr|o'f^(T) < oo and u(|x|) dx < oo. We are 
interested in the large-iV behaviour of the transformed path measure, 

e-HN,p-KN,p ( 1 , 3 ) 

Here f3 G (0, oo) is a finite time horizon which we will keep fixed in this paper. The trap part effectively 
keeps the motions in a bounded region of the space Through the pair interaction the i- 

th Brownian motion interacts with the mean of the whole path of the j-th motion, taken over all 
times before (5. Hence, the interaction is not a particle interaction, but a path interaction. We are 
most interested in the case lim^io^(^) = oo, where the pair-interaction repels all the motions from 
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each other (more precisely, their paths). In order to keep the notation simpler, we abstained from 
normalising the path measure in (USD. 

The model in Q was introduced and studied in [ARKOdj : see Section o for results from that 
paper and a discussion of the physical relevance of the model. In particular, in the limit f3 oo, 
followed by ^ oo, a certain variational formula appears that is called the Hartree formula in the 
literature. Therefore, we call the model in (dni) the Hartree model. 

When we take the limit as N ^ oo, we will not keep the pair-interaction function v fixed, but we 
replace it by the rescaled version = N'^~^v{N •). In other words, we replace with 

^ E ^ r r dsdt. (1.4) 

Note that N'^v{N •) is an approximation of the Dirac measure at zero times the integral of u o | - I, 
hence the double integral in Q is an approximation of the Brownian intersection local times at zero, 
an important object in the present paper. The Brownian intersection local times measure the time 
spent by two motions on the intersection of the their paths, see Section Id. II A natural sense can be 
given to this object only in dimensions d G {2,3}. Therefore, our analysis is naturally restricted to 
these dimensions. The main difficulty in the analysis of the model will stem from the Al-dependence 
and the high irregularity of the pair-interaction part. 

1.3 Results. 

We now formulate our results on the behaviour of the Hartree model in the limit as N —> oo, with /3 > 0 
fixed. First we introduce an important functional, which will play the role of a probabilistic energy 
functional. Define Ali(M'^) ^ [0, oo] as the Legendre-Fenchel transform of the map Cb(K^) 9/i—> 

^ logE[e-^o^-^^^'’^'^^] on the set Cb(M'^) of continuous bounded functions on where {Bs)s>o is one of 
the above Brownian motions. That is, 

Mfi)= sup ((/x,/)-ilogE[e/o"/(^«)d*]), /.gMi(]R‘'). (1.5) 

/eCb(R‘*)^ P ' 

Here Ali(M‘^) denotes the set of probability measures on Note that depends on the initial 
distribution of the Brownian motion. In Lemma 123 below we show that Jp is not identical to -|-oo. 
Alternate expressions for Jp are given in Lemma 12.31 below. Clearly, Jp is a lower semi continuous and 
convex functional on Ali(M'^), which we endow with the topology of weak convergence induced by test 
integrals against continuous bounded functions. However, Jp is not a quadratic form coming from any 
linear operator. We wrote {p., f) = J^dfix)p{dx) and use also the notation {f,g) = J^d f{x)g{x) dx 
for integrable functions f,g. If p possesses a Lebesgue density (jP' for some L^-normalised G L^, then 
we also write Jp^cfP) instead of Jp{p). In Lemma 12.21 below it turns out that Jp{p) = oo if p fails to 
have a Lebesgue density. 

In the language of the theory of large deviations, Jp is the rate function that governs a certain large 
deviation principle. (See Section [4.1 1 for the notion and some remarks on large deviation theory.) The 
object that satisfies this principle is the mean of the N normalised occupation measures. 


I ^ 

Pn,p = G N. 

i=l 

(1.6) 

p^^\dx) = ^J (5^(i)(dx)ds, i = l,...,N, 

(1.7) 
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is the normalised occupation measure of the i-th. motion, which is a random element of It 

measures the time spent by the i-th. Brownian motion in a given region. One can write the Hamiltonians 
in terms of the occupation measure as 

N 

= and KM,fS = l3 ( 1 - 8 ) 

i=l 

where we denote by V the integral operator with kernel v, i.e., V f{x) = vi\x — y\)f{y) dy and anal¬ 
ogously for measures Vfi{x) = J^d fJ-{dy) v{\x — y\). Hence, it is natural to expect that the asymptotic 
of the Hartree model can be expressed in terms of asymptotic properties of 

We are heading towards a formulation of our main result. Our precise assumptions on the trap 
potential, W, and on the pair-interaction functional, v, are the following. 

Assumption (W). W : ^ [0,oo] is continuous in {W < 00 } with limij^oo inf|x|>R = 00 . 

Furthermore, {W < 00 } is either equal to or is a hounded connected open set. 

Assumption (v). v: [0,oo) — > [0, 00 ] is measurable, Jjgdn(|x|)dx < 00 and n(|x|)^ dx < 00 . 


In order to avoid trivialities, we tacitly assume that the support of the initial distribution of the 
Brownian motions is contained in the set {W < 00 }. 

Now we formulate our main result. As we already indicated, the main role in the analysis of the 
Hartree model is played by the mean of the normalised occupation measures in m- 


Theorem 1.1 (Many-particle limit for the Hartree model). Assume that d G {2,3} and let W and v 
satisfy Assumptions (W) and (v), respectively. Introduce 


■= ^ [ v{\y\)dy 
Jw.d 


< 00 . 

1 ,,(i) 


(1.9) 


under the measure with density e with speed 


Fix P > 0. Then, as N ^ 00 , the mean = jy hp 0 / normalised occupation measures 

satisfies a large deviation principle on AAif} 

Nfi and rate function 

+ (W,(/>2) +4TTa{v) ||(/>||| if ^ exists, 

1 00 otherwise. 


i';Up) = 


( 1 . 10 ) 


The level sets {p, G Ali(M'^): < cj, c G M, are compact. 


We also write ii = ^. To be more explicit, the large deviation principle for means 

that 

• ] = “ weakly, ( 1 . 11 ) 

A—>00 iV p • 

where we identify Ali(M'^) with the unit sphere in L^(M'^) via the relation (/)^(x)dx = p{dx). The 
convergence in jnu is in the weak sense, i.e., the lower bound holds for open sets and the upper 
bound for closed sets (see Section HTTI for more details). Here we refer to the weak topology on Ali(M'^). 
See Section 13.11 for a heuristic explanation of the assertion of Theorem 11.11 

In Assumption (v) we require that n o | • | G L^(M'^). This is needed in our proof of the lower 
bound in (inii) only. We think that this assumption is technical only and could be relaxed if higher 
integrability properties of the elements of the level sets of were known. 

Interesting conclusions of Theorem II .11 are as follows. For a > 0, we introduce the variational 
formula 

XT{P)= inf (jp{P^) + {W,<l?)+4naM\i), 

0 eL 2 (Rd), || 0||2 = 1 V ^ / 


( 1 . 12 ) 
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which, for a = a{v), is the minimum of the rate function defined in (fTTUl) . Here are some facts 
about the minimiser in (EH. 

Lemma 1.2 (Analysis of Fix /3 > 0 and a > 0. 

(i) There exists a unique L'^ -normalised minimiser ( f >^ G L^(M^) n L^{W^) of the right hand side 

of (EH- 

(ii) For any neighbourhood N C n of (f^, 

Here ‘neighbourhood’ refers to any of the three following topologies: weakly in weakly in L\ 
and weakly in the sense of probability measures, if 4> is identified with the measure (j){x)'^dx. 

Now we can state some conclusions about the large-behaviour of the total expectation of the 
exponential Hamiltonian and about a kind of Law of Large Numbers. The proof is simple and omitted. 

Corollary 1.3. Let the assumptions of Theorem M . 1\ be satisfied. Then the following holds. 

(i) 

(ii) As N —>■ oo, piN,i 3 converges in distribution under the measure with density e 

towards the measure (^*(x)^dx, where 0* G L^(M'^) is the unique minimiser in (ITTl) with 
a = a{v) as defined in (EH)- 


1.4 Relation with quantum statistical mechanics. 

In this section we briefly explain the relation between the Hartree model in EHH and quantum statis¬ 
tical mechanics. 


An N-particle quantum system is described by the A^-particle Hamilton operator 
N 


Hm = 


Y^[-A‘f + W{x,) 


x = {xi,...,xn) G 


t,dN 


(1.14) 


2=1 


+ v{\Xi-Xj\ 

1<2<_7<7V 

For (d > 0 fixed, the trace of the Boltzmann factor, is given, via the Feynman-Kac formula, 

by a Brownian model similar to the one in El, where the Brownian motions are conditioned to 
terminate at their starting points (Brownian bridges) and the initial measure is the Lebesgue measure 
[br^ . However, the interaction Hamiltonian is, instead of A1 a,/ 3 in dl 




(1.15) 


This is a partiele interaction involving only one time axis for all the motions, in contrast to the 
time-pair integration in El- Note that there is no Hamilton operator such that the total mass of 
the Hartree model is equal to the trace of the corresponding Boltzmann factor. The Hartree model 
features the mutually repellent nature of the trace of a form which is more accessible to a 

rigorous stochastic analysis. The study of the large-behaviour of that trace, i.e., of the Brownian 
model with the interaction in (EH, is deferred to future work. 


For describing large systems of bosons at positive temperature, one has to consider the trace of the 
projection of Hat to the subspace of symmetric wave functions. The corresponding Brownian model 
is given in terms of Brownian bridges with symmetrised initial and terminal conditions. The effect of 
this symmetrisation on the large-limit is studied elsewhere |AKn6j . 
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Let us now comment on the physical relevance of the Hartree model. Its usefulness as a simplified 
model for the ground state of is well-known, see the physics monograph |DNn5l Ch. 12]. The precise 
relation is as follows. Recall the integral operator, V, with kernel v, and introduce the variational 
formula 


inf 


N 


2=1 


E (kivh])). 


(1.16) 


In |ABK04) we showed that there exist tuples (hi,..., h^) of minimisers, and we characterised them 
in terms of Euler-Lagrange equations and proved some regularity properties. Note that there is no 
convexity argument available, which leads us to the conjecture that the tuples of minimisers are not 
unique. Note also that 


Xn = iTf (^) H^h), where h = h\® ■ ■ ■ ® hjy. 


(1.17) 


Hence, one can conceive Xn' ns the ground product-state energy of H]\f, i.e., as the ground state energy 
of the restriction of to the set of N-fold product states. If ..., h)^') is a minimiser, we call 
:= (8) • • • <8) a ground product-state. One of the main results of |ABKn4j . Th. 1.7, states 

that, for any fixed G N, 

hm ^logE[e-^^./5-'^^-/5] = -xf. (1.18) 

p—>oo iV p 

The proof shows that the tuple of normalised occupation measures, (/u]g\ ..., (recall (ll.7j) l stands 
in a one-to-one relation with the minimiser tuples {hi,... of (11.161) . in the sense of a large 
deviation principle, analogously to (EU). This result illustrates the close connection between the 
zero-temperature Hartree model and the ground states of the Hamilton operator H^. One main 
novelty of the present paper is to demonstrate that also at positive temperature the Hartree model is 
a useful simplihcation of the trace of . 

It is instructive to compare the main result of the present paper. Theorem II .1 1 to the zero- 
temperature analogue of that result, which we derived in [ABKOdj and which served us as a main 
motivation for the present work. It turned out there that the Gross-Pitaevskii formula well ap¬ 
proximates the ground product-state energy x^^ in the large-N limit, provided that the interaction 
functional v is rescaled in the same manner as in dn. The Gross-Pitaevskii formula, derived inde¬ 
pendently by Gross and Pitaevskii in 1961 on the basis of the method initiated by Bogoliubov and 
Landau in the 1940ies, has a parameter a > 0 and is defined as follows. 

^ (llWll2 + (hk,())^)+47ra||())||^). (1.19) 

(/.eHl(Rd): ||0||2=1 ^ ^ 


There is a unique minimiser which is positive and continuously differentiable with Holder con¬ 

tinuous derivatives |LSSYn5| . 

The main result of the present paper. Theorem II .11 see (j1 .1 81) . is the positive-temperature analogue 
of the following result, which is |ABKn41 Th. 1.14]. 


Theorem 1.4 (Large-N asymptotic of Xn^)’ Let d G {2,3}. Assume that v satisfies Assumption (v). 
Replace v by vjsf{-) = N'^~^v{N ■). Let {h^p,... ,h^p) be any minimiser on the right hand side of 
(|1.16l) (with V replaced by vn{-) = N^~^v{N ■)). Define ■ Then we have 


hm 

N—^OO 


(®) ^ (GP) 

An 


and 


^N 


(Ci) 


where a{v) is the integral introduced in CH). The convergence of cfA is in the weak sense and 

weakly for the probability measures (/)^(x)da; towards the measure ((/>^°|jj)^(x) dx. 
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We remark that Lieb et al. proved the analogous result for the ground state of Hn, see Section ll^ in 
the appendix. Note that in d = 2 the scaling of v differs from the one used in Theorem II .41 Moreover, 
the parameter a{v) in Theorem II.41 is replaced by the scattering length of v in the result of Lieb et al. 
The integral a{v) is known as the first Born approximation of the scattering length of v [TaSSYD^] . 

We conjecture that in (|l-12l) converges to the Gross-Pitaevskii formula as —> oo. A proof 

of this is due to future work. 


2. Variational analysis 

In this section we derive some useful properties of the probabilistic energy functional Jp introduced 
in Q in Section EH and prove the existence and uniqueness of minimisers in the variational formula 
Xa’\P) introduced in (11.1211 in Section EH 

2.1 Some properties of Jp. 

First we show that Jg is not identically equal to +oo. 

Lemma 2.1. There is p, £ such that Jpip) < oo. 

Proof. Recall that Jp is the Legendre-Fenchel transform of the map Cb(]R^) 3 / logE[e^^bM/3>]^ 

where we recall that p/s is the normalised occupation measure of one of the Brownian motions. Recall 
the mean of the N normalised occupation measures from (|1.6I) . Now pick a continuous function 
—> [0,oo) satisfying limij^oo inf|a;|>R5(a^) = oo- Then we have, for any (7 > 0, by splitting the 
probability space into {{g,'PN,f3) ^ C} and its complement, 

—00 < logE[e“''^’^'^^] = limsup — logE[e“'^^^’^Y/3)j 

N^oo -V 

< max I - G, limsup^ log I. 

According to [DZ981 Th. 4.5.3(b)], the sequence {'PN,i 3 )NeN satisfies the upper bound in the large 
deviation principle for compact sets with rate function equal to Jj^. By Prohorov’s Theorem, the set 
{p G Ali(R'’*): {g,p) < G} is compact. Furthermore, note that the map p —{g,p) is upper semi 
continuous. Hence, the upper-bound part in Varadhan’s Lemma, |DZ981 Lemma 4.3.6], implies that 

limsup ^ log E [e“^<^’^Y/3) < - inf, ^ ({g,h) + Mh'))- 

Picking G large enough, we find that 

00 > inf, , ((5 ,r) + «//3(r))- 

This implies that Jg is not identically equal to 00 . □ 

Now we show that Jg is infinite in any probability measure in R'’* that fails to have a Lebesgue 
density. 

Lemma 2 . 2 . If p £ Ali(R'’*) is not absolutely continuous, then Jfsip) = 00 . 

Proof. We write A for the Lebesgue measure on R'^. Pick p £ Ali(R'^) that is not absolutely 
continuous. Then there is a Borel set A C R'^ such that A (A) = 0 and p{A) > 0. Let M > 0. We show 
that Jf^{p) > M. Pick K = and g = ^p{A). We may assume that 2 < {Q£)e>o be 

an increasing family of open subsets of R'’* such that A C Qe and XiQs) < e for any e > 0. Let pp is 
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the (random) normalised occupation measure of a Brownian motion {Bs)s>o- Pick e > 0 and <5 > 0 so 
small that 

where Us{Qe) is the ^-neighbourhood of Q^. This is possible since 


limsup limsup/i;3(?75((5e)) = limsup 
ejO <5t0 ejO 


(n<?, 


£>0 


= 0 a.s. 


Now we pick a function / G Cb(I^^) such that 0 < f < K , supp (/) C Us{Qe) and /|q^ = K. Then 
(/ih) ^ /qj/(®)> Kfi{A). Furthermore, 


E 


= E 


< I?} +E > v} 


+ e^^F{^,0{Us{Qe))>v) 

< qI^vK g/iRg-ZiRii-*?) < < QpKlt,{A)_ 


( 2 . 1 ) 


Hence, 


Jjiilj) > (/,A^) - ^logE 




> Kix{A) - K-n{A) = M. 


□ 


In general, the supremum in the definition of Jjj is not attained. It is of interest to replace 
the function class Cb(IK‘^) in by some class of better behaved functions. In particular, one would 
like to use only functions / that are extremely negative far out. This is of course possible only if cj)^ 
decays sufficiently fast at infinity. We write for the normalised occupation measure of the Brownian 
motion {B^p)s>o in the following. By m we denote the initial distribution of the Brownian motions. 

Lemma 2.3 (Alternate expression for J^). Let W : —> [0,oo] be continuous in {W < oo}, which 
is supposed to contain supp(m) and to be either equal to or compact. Fix cj) G L^(]R‘^) satisfying 
{W, 4>‘^) < oo. Then 

sup ((-lT + /i,<^2^-4logE[e^<-'^+'^’^^>]). (2.2) 

/leCbfK'*) ^ P ^ 


Proof. Let Jw,p{4>‘^) denote the right-hand side of (12.21) . We first prove ‘>’ in (12.21) . Let h G Cb(lK'^)- 
We may assume that h < 0 (otherwise we add a suitable constant to h). Then fn = (—W-l-h) V (—ii) is 
a bounded continuous function for any i? > 0 with (—IF+/i, fp‘) < {Jr, fp‘). Furthermore, /r J, —IF +h 
as i? —> oo, hence the monotonous convergence theorem yields that 

lim ^ Elg/Jf-VP+h./x;,)! _ 

R^oo L J L J 

This shows that Jp{4>^) > Jw,(5{4>^) holds. Note that we did not need here that fW^cjP') < oo. 

Now we prove ‘<’ in (12.21) . Let / G Cb(M'^) be given. For i? > 0, consider Lr = {f + IF) A R, then 
Lr G Cb(l^^) with hi? I / + IF. Since (IF, < oo, we have 

liminf(-IF + /ir,(/>^) > {f,c/)^). 

R^oo 

Furthermore, by the monotonous convergence theorem, 

lim . 

R^oo J L J 

This implies that ‘<’ holds in dH. □ 

Let us draw a conclusion for compactly supported functions (j). For a measurable set A C we 
denote by Cb(^) Ike set of continuous bounded functions A —> M. 
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Corollary 2.4. Fix cj) G satisfying ||<^||2 = 1- If the support of cj) is compact, connected and 

contains supp(m), then 

sup - ilogE[e^<^’^^>l|3upp(^^)csupp(<«}])- (2.3) 

/eCb(supp(</i)) ^ P ^ 

Proof. We pick W = oo]lsupp((/i) iii Lemma ITHl and see that, on the right hand side of (12.21) . we may 
insert the indicator on {supp(/i_( 3 ) C supp(()))} in the expectation and can drop W in the exponent. 
Hence, both this expectation and the first term, {—W + f,(j>^), do not depend on the values of / 
outside supp((/>). □ 

The next lemma shows the interplay between the arguments for the functional Jp and the fixed 
initial distribution of the Brownian motions. 

Lemma 2.5. Let 4> G satisfying \\4>\\2 = 1- //dist(supp((^), supp(m)) > 0, then = +oo. 


Proof. Let S be an open neighbourhood of supp(i?!)) with 6 = dist (S', supp(m)) > 0. Pick K > 0 
and a continuous bounded function /: —> [0,K] with supp(/) C S and /Isupp(fli) = Then 

(/j = K- Then we have 

E[e/3(/-M/3)] = E[e^T,/^/3>]l{/i;3(5) < 1 - > 1 - 

' I 


Hence, 


sup \Bt\ >6 

'0<l<A-i/2 


< e' 


I3K 






> (/,/> - ilogE[e^(/«)] > _ 1 i„g 


/? 




Letting iL —> oo shows that Jg{4>‘^) = +oo. 


(2.4) 


□ 


2.2 Analysis of Xa^iP)- 

Proof of Lemma ll.2L The uniqueness of the minimiser follows from the convexity of the functionals 
Jp and {W, •), together with the strict convexity of fP e-> \\4>\\\- 

Let ((^n)nGN be an approximative sequence of minimisers for the formula in (11.1211 . i.e., cfn £ L^(M'^), 
|||?^n ||2 = 1 for any n G N and 

+ (bL,0n) +4vra UnWtj = Xa\P)- 

In particular, the sequences (T/ 3 (<?i^))n) ((kP,(/)^))„ and (||(/>n|| 4 )n are bounded. Since W explodes at 
infinity by Assumption (W), the sequence of probability measures (^^(rE) da;)„£is} is tight. According 
to Prohorov’s Theorem, there is a probability measure p, on such that (jf^ (x) dx converges weakly 
towards as n —> oo, along a suitable subsequence. Since the sequence is bounded in L^(R'’*), the 

Banach-Alaoglu Theorem implies that we may assume that, along the same sequence, (<?i>^)n converges 
weakly in L^(R'’*) towards some (jP G L^(R'’*). Since Jp is weakly lower semi continuous (in the sense 
of probability measures), we get 

liminf J^(())2) > Jpip), 

n^oo 

and with LemmaE^Jwe conclude that p{dx) = dx for some function (jp G L^(R'’*) with ||<?i ||2 = 1- 
By the weak convergence in L^(R'’*), combined with the weak convergence in the sense of probability 
measures, for any continuous bounded function V’ with compact support we have 
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Hence, we get cfP' = (jP' a.e.. As Jp is weakly lower semi continuous (in the sense of probability 
measures), || • ||4 is L^-weakly lower semi continuous, and (j) {W^cfP) is lower semi continuous, we 
have that 

+ {W,4P) + 'i'Kaiv) \\(l)\\\ < Iminf (jp{(lPn) + + 4vra(n) \\(t>n\\fj 

This shows that the limiting point (j) is the minimiser in (11.121) . This ends the proof of Lemma ll.2r il. 
Furthermore, the proof also shows that the minimising sequence converges, along some subsequence, 
towards the unique minimiser in all the three weak senses: in L^, and weakly as a probability 
measures. This implies Lemma nn ii). 

□ 


3. Large-behaviour: Proof of Theorem 11.11 


In this section we prove Theorem 11.11 We shall proceed according to the well-known Gartner- 
Ellis theorem, which relates logarithmic asymptotic of probabilities to the ones of expectations of 
exponential integrals. Therefore, we have to establish the existence of the logarithmic moment gen¬ 
erating function of Tiv ,/3 under the measure with density e 
Theorem II .Il ls the following. 


AN) 


M/3. The main step in the proof of 


Proposition 3.1 (Asymptotic for the cumulant generating function). For any f £ Cb(]R'^) the cumu- 
lant generating function exists, i.e., 


where 


X" 


lini log E 
A—>oo N(3 


*(/) = ~ inf 

P </.eL2(Kd): I 


(JV) 


' ^ N,0q^ if Wn,0) 




+ {W - f, fP) + 47ra{v)\\4>\\l^ 


(3.1) 


(3.2) 


Indeed, Theorem II.II follows from Proposition 13.11 as follows. 

Proof of Theorem urn We are going to use the Gartner-Ellis Theorem, see |DZ981 Gor. 4.6.14]. 
Eor this, we only have to show that the sequence of the ]In,p is exponentially tight under measure 


-H 

with density e n ,0 and that the map / is Gateaux differentiable. 

The proof of exponential tightness is easily done using our assumption that lim/j^oo inf|a,| >rW{x) = 
oo in combination with the theorems by Prohorov and Portmanteau; we omit the details. 

Fix / £ Cb(lK'^)- The proof of Lemma EH shows that the infimum in the formula of the right hand 
side of (El is attained. Let (jPj: £ L^(R'^) be the minimiser for the right hand side of ra, and, for 
some g £ Cb(lR^)) let fPfj^^tg ^ T^(]R'^),t > 0, be corresponding minimiser for f + tg instead of /. We 
obtain 

+ >{9.'/L„). (3.3) 


X' 


Since {(lPf^tg)t>o Is easily seen to be convergent weakly (in the sense of probability measures) towards 
the minimiser (jpj of the right hand side of O, it is clear that the right hand side of (El con¬ 
verges towards {g, Analogously, one shows the complementary bound. This implies the Gateaux- 
differentiability of with 


A 

dg 


x^®H/) = ( 5 , <(>/)• 


□ 
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In Section 123 we give a heuristic explanation of (EH) and introduce the Brownian intersection local 
times, an important object in our proof. In Sections III.21 andrespectively, we prove the upper and 
the lower bound in (EU). 


3.1 Heuristics and Brownian intersection local times 

In this section, we give a heuristic explanation of the assertion of Proposition I3.1L We rewrite the 
two Hamiltonians in terms of functionals of the mean defined in (USD and use a well-known large 
deviation principle for particular, we introduce an object that will play an important role in 

the proofs, the Brownian intersection local times. For the definition and the most important facts on 
large deviation theory used, see the Appendix or consult |DZ98j . 

Rewriting the first Hamiltonian in terms of /Utv,/? is an easy task and can be done for any fixed N: 

1 ^ 

Hn,p = NP / W{x)-J2^^0idx) = Np{W,JI^^p). (3.4) 


Now we rewrite the second Hamiltonian, which will need Brownian intersection local times and an 
approximation for large N. Let us first introduce the intersection local times, see |(IHR.84j . For the 
following, we have to restrict to the case d G {2, 3}. 

Fix 1 < i < j < N and consider the process Ri'i — the so-called confluent Brownian motion 
of Ri'i and —B^^K This two-parameter process possesses a local time process, i.e., there is a random 
process {L^^’^\x))^^^d such that, for any bounded and measurable function /: —> M, 


/ fix)L^g'^\x)dx 

jRd 

Hence, we may rewrite 


= ^ [ ds / dt f {B^fl - Bl^'>) = f j //J;’^'\dx)^J;’"'\dy)/(x-y). 
P Jo Jo Jw^ Jw^ 

as follows: 

^ [ vizN)L^^’^\z)dz 

= NP f v{x)^ ^ dx. 


(3.5) 


(3.6) 


It is known |OHR,84l Th. 1] that (®))xgR‘^ may be chosen continuously in the space variable. Fur¬ 
thermore, the random variable Lp’^\0) = \\m.x^o L^lP^\x) is equal to the normalised total intersection 
local time of the two motions and B^^^ up to time fl. Formally, 




/ / 

Ja Jo 


^ , l{Bi'^ G dx} 


ds 


dx 


I 


^ t{Bf G dx} 
dx 


Using the continuity of we approximate 



/ig\dx) ^if{dx) 
dx dx 


(3.7) 


N,i3 


Npd7ra{v) L^id) 


N pdTTa{v) 


1<2<J<A^ 


N N 


2=1 


2 = 1 


N P'iTra{v) 


dpjv,/t 2 

dx 2 


where we conceive as densities, like in (I3.7jl . 

The main ingredient is now that (77 a,/3)a^gn satisfies a large deviation principle on AIi(M'^) with 
speed Nfl and rate function Jp. This fact directly follows from Cramer’s Theorem, together with the 
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exponential tightness of the sequence (^jv^^)ArgN- Hence, using Varadhan’s Lemma and ignoring the 
missing continuity of the map // i—> ||^|| 2 ) this heuristic explanation is finished by 


E 




(JV) 


E 


exp<^ -Nl3 


(VL -/,/xjv,/3> - 4vra(u) 




dx 

2 J j. 




Here we substituted (/>^(x) dx = /u(dx) and noticed that, according to Lemmawe may restrict the 
infimum over probability measures to the set of their Lebesgue densities (jy^. 


3.2 Proof of the upper bound in 1|3.1|1 in Proposition 13.11 

In this section we prove the upper bound in 13.II in Proposition I3.1L Our proof goes along the lines 
of the argument sketched in Section 13.1 1 However, in order to arrive at a setting in which we may 
apply Cramer’s Theorem and Varadhan’s Lemma, we will have to prepare with a number of technical 
steps. More precisely, we will have to estimate the interaction term from below in terms of 

a smoothed version of the intersection local times. This version will turn out to be a bounded and 
continuous functional of the mean of the normalised occupation times measures, which are the 

central object of the analysis. 

Our strategy is as follows. First, we distinguish those events on which, for at least ((1 — r])N)‘^ 
pairs (i, j) of indices, the intersection local times L^“’^'(x) for |x| < 2e are sufficiently close to 
and its complement. More precisely, we will have |L^*’^'(x) — ^ for these {i,j) and x. Here 

rj, e are positive parameters which will eventually be sent to zero. (The complement of the event 
considered will turn out to be small by the continuity of the intersection local times in zero.) The 
replacement of by will require a space cutting argument, i.e., we will restrict the 

interaction from to the cube Qji = [—R, i?]^ for some R > 0 which will eventually be sent to 
infinity. Our second main step is to replace by the smoothed version * Ks{0), where 

Ks is a smooth approximation of the Dirac function (5o as e | 0. The smoothed intersection local times 
can easily be written as a continuous bounded functional of the mean of the normalised occupation 
times measures, Hence, Cramer’s Theorem and Varadhan’s Lemma become applicable, and 

we arrive at an upper bound for the large-V rate in terms of an explicit variational formula, which 
depends on the parameters. Finally, we send the parameters to zero and infinity, respectively. 

Let us come to the details. Introduce the following random set of index pairs, 

= = sup |L(-\0)-4 -\x)| <d, (3.8) 

|a:|<2e 

where = {{i, i): i G {1,..., A^}} denotes the diagonal in {1,..., A^}^. Fix ry > 0 and consider the 
event 

An = ANie, e, r?) = {3/ C {1,..., V} : I x / \ C Dn, \I\ > (1 - r?)A^} (3.9) 

In words, on An, there is a quite large set I of indices such that all pairs {i,j) of distinct indices in I 
satisfy — L^’^'(x)| < ^ for all |x| < 2e. 

First we show that the contribution coming from the complement A^ vanishes for small e: 


Lemma 3.2 (H^ is negligible). For any ^ > 0 and any rj G (0, ^), 


lim lim sup — log E 

£—^0 J\J'—IV 


-H. 






— OO. 


(3.10) 


Proof. 


Since HN^fi, bounded from below, it suffices to show that 


lirn lim sup 4 log P (H^ (e, C, r/)) 
e^O N^oo JV 


= —OO. 


(3.11) 
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Note that 

A% = A%ie,^,v) = {yic{l,...,N}-. \I\>{l-rj)N^{IxI)\AN<^DN}. (3.12) 

First we show that, on Aj^, there are pairwise different integers *i, ji, *25 j2) • • • j *[r 7 Ar/ 3 j > Jl» 7 A/ 3 J iii 
{1,... ,N} such that G for alH = 1,..., [r]N/3\. 

We construct the indices inductively. Consider Ii = {1,... , [(1 — 17 )^']}, then 

/i X Ii \ Atv ^ Dn, 

i.e., there is a pair (ii, ji) G {h Ii\ ^n) C Also the set 

h = ih \ {kJi}) U {(1 - v)N + 1, (1 - v)N + 2} 

has no less than (1 — r])N elements. Hence there is a pair G {I 2 x I 2 \ ^n) C 

Clearly ji, 12 ,^ 2 } = 4. In this way, we can proceed altogether at least times. 

This procedure constructs the indices ji, ^2j • • •, *[r;7v/3j > J[r?7V/3j pairwise different such that 

- , (j[r?Af/3j )iL»?A/3j) G -C’n- 

Now we prove that (TOT]) holds. We abbreviate ‘p.d.’ for ‘pairwise disjoint' in the following. For 
notational convenience, we drop the brackets [-J. Because of the preceding, we have 

F{A%iE,^,rj)) < = 1’ • • • sup ^ 

■ ■ ^ ATX ^ ^ kl<2e / 



(3.13) 

Since the process iLp’^\x))^^^d may be chosen continuously in the space variable |CHR.841 Th. 1], we 
have 

limPf sup \Lk‘^\x) — = 0. (3.14) 

<^^0 ^|a;|<2£ p p J 

This, together with (EH, concludes the proof. 

□ 

Now we estimate on the event Ajq. For i? > 0, we recall that Qr = [—R,R]‘^ and introduce 

aR{v) =-^ [ v{\x\)dx. (3.15) 

8 ^ Jqr 

Let k: —> [0, 00 ) be a smooth function with support in [—1, l]'^ and f k(x) dx = 1. For e > 0, we 

define Ki;{x) = e~^K{x/e). Then Kg is an approximation of do as e | 0, and we have supp <G Qe and 
Ks(x) dx = 1 for any e > 0. 

Lemma 3.3 (Estimating on AN{e,^,r])). Fix e,^,r] > 0. Then, for any R> 0 and any A G N 
satisfying N > R/{2e), on the event Aiq{e,f^,ri), 

-A'^^ < -d'KaR{v)f3\I\ (1 - ri)\\Ui,(i * i^e\\\ + 8TTf,aR{v) + l°° , (3.16) 

where the random subset I of {!,...,A} in (13.91) is chosen minimally with |/| > (1 — f])N and 
(/ X I \ Ajv) C Dn, and 

' ' i&i 

denotes the mean of the corresponding normalised occupation measures. 


(3.17) 
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Proof. First, we write the interaction terms for the scaled pair potential vn as integrals against 
the intersection local times of two Brownian motions at spatial points x/N. As the pair interaction 
V is positive we get easily an upper bound, when we restrict the integrations to the box Qr. On this 
compact set we use the continuity to get rid of the dependence of the spatial variables in the integrals. 
We restrict the summation over all pairs of indices to the random set Djy = L)jv(e,0- 

Recall JSSl). On the event Ajv(e,C,i?), we may estimate, for all N > R/{2e), 




1<2<ji<7V 




da; 


X 


< —dvr 


OiR{v) 

iV2 




{id)&DN {i,j)eDN 

<2^ Ly>{{))+ 87raR(u) 


dx 


(3.18) 


Ar2 




< yy Ly\0) +87r^aR{v). 


iV2 


{i,j)^Disr 


For any (i, j) £ D^, we now replace the intersection local time at zero, L^’^^(O), with the smoothed 
version * Ke(0). The replacement error is estimated by 


Ly \0) - * Ks* Ke(0) 


(i,i) 


< 


dxKe * Keix) Ly\0) - Ly\x) < (i, j) G Dm- (3.19) 


-(*.j) 


' Q2e 


Hence, we may continue the estimation in (mi by 

Nf5^^N,i3 - 


'yy ly^ *Ke* Ks{0) + 8TTCaR{v). 


(3.20) 


{iJ)€DN 

Recall the defining property on the intersection local time in dsn. Hence, we can write the smoothed 
version in terms of a convolution of the normalised occupation measures as follows: 


Ly'' * * ^^(o) = 


( 0 ) = / / / ^J-^y(^x)yy{dy)dzK^{y-x- z)ks{z) 

jRd jRd jRd 

= [ [ [ yf{dx)yy{dy)dwKs{w-y)K,{w-x) 

jRd jRd jRd 

= * He)- 


(3.21) 


To proceed with our estimates from (mi, we add now the self-intersection terms, i.e., the diagonal 
terms where i = j. The additional terms are bounded by the L°°-norm of as seen from 


{^3 * ^ 6,^3 * He) = [ [ f dx yydy) yy{dz) Keix - y)Hsix - z) 

jRd- jRd jRd- 


< \ \k. 


— 11 11 oo 


Hence, we obtain from (TOl that 


-^K[y < -471 


NP 


N,f3 - 


ORiv) 


[ [ dxyydy)K,{x-y) = \\KPU 

jRd jRd ^ 

aR{v)\\Ke 


(3.22) 


yy {^0 * He, yy * ny 8Tr^aR{v) -h 477- 


(i j)G£*ivUAjv 


N 


(3.23) 


Recall that we work on the event Atv defined in (mi). On this event, Dn U Aat contains a subset 
of the form I x I with I C {1 ,... ,N} and |I| > (1 — r])N. Let such a random set be chosen, for 
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definiteness we choose it minimally. We continue the estimation of (non with 


1 


aR{v)\I\ 


aR{v)\\Ke 


^ * ^s) + 8.CaR{v) + 

and from this the assertion follows. 

Using Lemma lOl on the left hand side of m, we obtain the following bound. 
Corollary 3.4. Fix e,^,rj > 0. Then, for any R > 0, as N ^ oo, 


E 


exp|-/3[(l - r])N\ (^{W - f,lIi(i-rj)Ni,f3) + 47raij(u)(l - v)\\Tl{i-n)N\,0 * Z\T) } 


< Qo{N)^Nl'nf}\\f\\oo+8ihaRiv)+Cv] 


X E 


where C > ^ is an absolute eonstant. 


(3.24) 

□ 


(3.25) 


Proof. Since the trapping potential W is nonnegative, we easily estimate 

N 

i=l iel 


(3.26) 


Furthermore, it is easy to see that 

Nf3{f,-p^^^)<P\I\{f,Tlj^f,)+N7jP\\f\U (3.27) 

Using these two estimates and (ixm in the expectation on the left hand side of (EU), we obtain 


E 


^N,SQ^{fWN,l3)^ 


An{£,^,v) 


< E 


exp|-/3|/| (^{W - f,ni^p) + iTraR{v){l - v)\\Ti,f} * ^ell^) }^^iv(£,L 


v) 


,,ivc„,e.fl+o(N) 


(3.28) 


where 


Cr,,^,R = vP\\f\ loo + S-Kij3aR{v). (3.29) 

Now we sum over all possible values of the random set I and note that the distribution of is equal 
to the one of Ri p = j X]!=i Z ■ Hence we get 

l.h.s. of E2H1) < e^(^r,x,R+oi.N) Y 

l=\(l-ri)m Lc{l,...,Ar}: \L\=l 

E l/=Lexp| +47raij(u)(l 

N 

< gNC„,5,fl+o(A) Y 


X E 


N 
. I 

«=[(1-7?)AJ 

exp I - ip(^{W - + A-KaR{v){l - ri)\\jlip * k, 


||2 
'e|l2 


)} 


(3.30) 

It is clear that there is a C > 0 such that, for any ly > 0 and any I E {[(1 — r])N~\,..., N} we 
can estimate ('^) < . In the exponent, estimate I > [(1 — i?)iVj where I is multiplied with a 

nonnegative factor, respectively estimate l{f-,Tii^p) < [(1 — ri)N\ + ??iV||/||oo- The sum on I is 

estimated against N, which is absorbed in the term e°('^^. □ 

Now we use arguments from large deviation theory to identify the large-rate of the right hand 
side of (E2S1): 
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Lemma 3.5 (Large deviation rate). For any > 0 and any R> 0, 


1 


limsup-—logE 


< r](3\\f Woo + 8TT^(5aR{v) + Crj - {I - r])xi%^r^{f), 


(3.31) 


where 


x?i,(/) = ^ ^ 4„<,„(.,)(1 - ,)||02 . „.||i}, (3,32) 


Proof. On the right hand side of (I3.25|) . we will apply Cramer’s Theorem jDZ981 Th. 6.1.3] in 
combination with Varadhan’s Lemma |DZ981 Lemma 4.3.6] for the mean of normalised occupation 
measures. Let us explain what these theorems say and how we apply them. See the Appendix, 
Section EH for a brief account on large deviation theory. 

Note that Ali(M'’*) is a closed convex subset of the space of all finite signed measures, A1(M'’*). The 
duality relation 

(/, p) G Cb(M'^) X M(M'^) ^ [ f{x) n{dx) (3.33) 

jRd. 

determines a representation of A1(M'’*)*, the topological dual of A4(M^), as Cb(M'^). Clearly, the 
topology inherited by Ali(M'’*) from A1(M'’*) is the weak topology, which is induced by the test integrals 
against all bounded continuous functions. Ali(M‘’*) is a Polish space with the Levy metric |DSnij . 
Thus, all the assumptions of [DZ981 Th. 6.1.3] are satisfied. Hence, {JIi^i 3 )i£N satisfies a weak large 
deviation principle. The rate function is equal to the Legendre-Fenchel transform of the logarithmic 
moment generating function of that is, it is the function Jp defined in (11.51) . 

For any M > 0, the set 

Fm = {fie Mi(M'^); {W,fi) < M} 

is compact, as is easily derived with the help of Prohorov’s Theorem, using that lim|3,|^oo W(x) = oo 
and the lower semi continuity of the map fi e-> (IT, fj,) (see Assumption (W)). 

We abbreviate I = [{1 — r])N\ for a while. On the right hand side of H3.25I) . we insert + 

On the event {]ii,p e F^j}, we estimate the trap part, {W,]Ipp), from below against M and 
use that the remaining terms in the exponential are bounded from below. Hence, 


E 


exp^-lp(^{W - f, fipp) + 4TTaR{v){l - r])\\fii^p * 

< + E exp|-//3(^(lT - f,JIi,p) + 47 raij(u)(l - v)\\'Fi,i 3 * 


(3.34) 


The functional fi i—> {W — f, fi) + d7raR{v){l — rf)\\fa*K £\\2 is lower semi continuous, as is easily seen from 
Fatou’s lemma, and bounded from below. Furthermore, according to |DZ981 Th. 4.5.3], the family 
i'Fi,i 3 )i£N satisfies the upper bound in the weak large deviation principle with rate function Jp. Hence, 
we may apply the upper-bound part of Varadhan’s lemma |DZ981 Lemma 4.3.6] to the right-hand side 
of H3.34I) . to obtain, if M is sufficiently large. 


lim sup — log E 
1^00 


exp 


|-Z/3(^(1T - f,fipp) +A-KaR{v){l - rf)\\fiip * K^Ha)} 

v) + {W - /, v) + d'naR{v){l - ??)||i^ * • 


< — inf 
i'GFm 

Lemma 12.21 implies that the infimum on the right hand side is equal to X^^Rrj 
Summarising the contributions from Ari and A^, we obtain the following. 


if)- 


(3.35) 


□ 
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Corollary 3.6. 

lim sup 


N^OO 


NP 


logE 


(N) 




< — lim inf lim inf lim inf XeR nif)^ 
?jJ.O R—>oo ej,0 ’ ’ ' 


(3.36) 


where X^^Rnif) defined in (j3.32B . 


Proof. A combination of Lemmas EJ and 13.51 gives that the left hand side of is not smaller 

than 

vP\\f\\cx> + ^PaR{v) + Cri - (1 - 7?) liminfxf’ (/), 

eiO ’ 

for any rj,R,^ > 0. Letting r/, ^ | 0 and R —> oo, we arrive at the assertion. □ 

Now we identify the right hand side of (E2SI): 

Lemma 3.7 (Approximating the variational formula). 

lim inf lim inf lim inf x^®p(/) > x^®^(/), (3.37) 

»?i0 R^oo ejO ’ ’ ' 

where X^®H/) defined in H3.ll) . 


Proof. 


We hrst fix r/ > 0 and R > and prove that 


lim inf xS,^(/) > Xr^, 


(3.38) 


with 

+ {W - /,(/>^) + 47raij(i;)(l -r/) 110114 ). (3.39) 

Let (0e)£>o be an approximate minimising sequence for the right hand side of (j3.39|) . i.e., 

^(^/ 3 ( 0 £) + {W - /, 02 ) + 47 raR(u)(l - 77 )] 10 ^ * = li^inf x^®^,,?(/)• 

In particular, (Jg(0^))£>o, ((IP — /, pe))£>o and (||0g * Ke\\‘^)e>o are bounded. As e | 0, along suitable 
subsequences, the probability measures 77 e(dx) = 0 g(x)dx converge weakly to a probability measure 
^{dx). Certainly, we may assume that X^f^Rnif) bounded as e | 0, and therefore also Jpicj)^) is. 
The lower semi continuity of Jg with respect to the weak topology of probability measures gives that 
< 7/3 (/i) < lim inf£\o < 7/3 (lie) < 00 . From Lemma ^I?]\ we get the existence of a density for the measure 
17 , i.e., li(dx) = 0^(x)dx for some cjfi G L^(M'^) satisfying || 0||2 = 1. Since e ^ ||0e * i^elb is bounded, 
there is a function 0 G L^(M'^) such that, along suitable subsequences, 0 ^ * converges weakly in 
L^(M'^) to and II 0 II 4 < liminf^io ||0e * Thus, for every g G Cc(M'^) we get 

I { 9 , 0^ - 0^) I < I {g, 0 ^) - ( 5 , 0^) I + I ( 5 , 0£ ) - ( 5 , 0£ * Ke) I + I ( 5 , 0£ * Ke - 0^) I • (3.40) 


Now, the first term on the right hand side of (lOnl) vanishes in the limit e —> 0 because g is bounded 
and continuous and the convergence follows from the weak convergence of the probability measures. 
The second term on the right hand side of (TTini) is estimated as 

\{9AI) - * Ke)\ = \{9- 9* K£,0£)| <\\g-g* Kelloo ^ 0 for e ^ 0, (3.41) 

since g G Cc(M'^). The last term on the right hand side of (mni) vanishes in the limit e —> 0 because g G 
L^(]R'^). Hence, {g,4>‘^) = (5,0^) and therefore 0 = 0 almost everywhere. Clearly, liminfe|o(kF,0^) > 
{W^cfP). Altogether, (I3.38j) follows. 

We now finish the proof of the lemma by showing that 

lim inf lim inf x^ ’ (/) > x'®’ (/) • (3.42) 

riiO R^oo ’ ' 

Since the map x fda^l) is assumed integrable, it is clear that, as 72 —> oo, aR{v) converges towards 
q;( 7;) dehned in (Hini). Hence, the proof of is an easy task, and we omit it. □ 
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3.3 Proof of the lower bound of (|3.1 jl in Proposition 13.11 

Now we turn to the proof of the lower bound in (EU). We were not able to produce a proof along the 
lines of usual large deviation arguments including Cramer’s Theorem and Varadhan’s Lemma, since 
we did not find any way to overcome the technical difficulties stemming from the singularity of the 
pair interaction term. Instead, we write the expectation on the left-hand side of (EU) as N iterated 
expectations over the N Brownian motions and use a lower estimate that is directly implied by the 
definition of the rate function, J^, more precisely, of some variant to be introduced below. In order 
to explain this idea, fix (/> G L^(]R'’*) and consider the random potential 

qr-=- {N - j)VN * j = (3.43) 

i=l 


Here we used the notation V * v{\x — y|) fj,{dy) (analogously with (j)^, conceived as a finite 

measure) and an analogous notation for V]\f with v replaced by ujv(') = ■ N). We rewrite the 

left-hand side of (ED) as follows. We write for the expectation with respect to the Lth Brownian 
motion. Recalling the definition of in (11.41) and noting that qn = — we see that 

we have 


E 




= E^^^ (g) • • • 




X E^^^ —W+fp 


{N)s 


(3.44) 


The main idea in our proof of the lower bound is that the definition of directly implies the estimate 


^ h E Cb(M'^),,/> E Uh = L (3.45) 

If the random potential h = — W + f were in Cb(E'^) almost surely, then (1X151) instantly implied 

an estimate for the last term on the right hand side of (IXIH) . and we could choose (j) arbitrary and 
optimise later on cj). 

However, the random potential h = qi\f — W + f does not have sufficient regularity for applying 
(1X151) directly. But note that qj lies in L^(M'^) almost surely, as is easily derived from the assumption 
that Jjgd t'(|aj|)^ dx < oo. In order to make (1X151) applicable for functions h of the form q — W + f 
with q E L^(M'^), we have to establish a lower bound for in terms of a supremum over this class of 
potentials. This is the content of the following lemma. 

Lemma 3.8. For any (j) £ L^(M'’*) IH L^(M‘^) satisfying \ \4>\\2 = 1; 

sup ((-IT + h,(^2^ - ilogE[e^<-^+^’'^/3>]). (3.46) 

hGL2(R'^): /i<0 Y 2 


Proof. As a first step, we show that 

sup ((/,(/)^) - ilogE[e^<^’'"^>]). (3.47) 

/GLRR'^): /<0 ^ y ' 

Let / E L^(M‘’*) be given satisfying / < 0. We approximate / by continuous bounded functions in a 
standard way as follows. Let k: M'’* —> [0,oo) smooth with n(x) dx = 1 and supp k C [—1,1]'^. Fix 
e > 0 and consider fe = f * where ^^(x) = e~^n{x/e) for x E M'’*. Then f^- is continuous. Using 
Schwarz’ inequality and the fact that / E L^(]R'^), one sees that is bounded. Furthermore, 

liminf(/£, = hminf(/, (j)^ * k^) = (/, (3.48) 


since / E L^(R'’*) and (j)^ E L^(M‘’*). 
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Since fe—^f strongly in we can pick a subsequence Sn i 0 such that fe^ - 
everywhere. Since / < 0, it follows from the bounded convergence theorem that 

lim ETe-^o^ feniBs)dsi = jEle/o'' fiBs)dsi _ 

n—>oo 

From this, together with (ITlHll . (HTITI) follows. 

Now we prove by showing that 

1 

P 


f pointwise almost 
(3.49) 


> sup ((-kF + /i,(/>2) - ilogE[e^<-^+'^’^/3>]y 
h£L^{Rd.): h<0^ P ' 


sup ((/,,^2^-llogE[eW.M]) 

/GL2(Rd), j<o \ P 2 

This is similar to the proof of ‘>’ in O in Lemma El Let h G satisfy /i < 0, and consider 

fu := {—W + h)tQj^ for R > 0. Clearly, fji G with fji < 0, and /r | —W + h pointwise as 

i? —> oo. Therefore, 

liminf(/R,</>2) > {-W + h,(tP), 

R-^oo 

according to the monotonous convergence theorem. Furthermore, for the same reason, and since 
—W + his bounded from above, 

limsupE[e^<-^«’^/3>] < 

i?—>C» 

This implies the statement and finishes the proof of (EH- □ 

We would like to remark that it is the assertion in that forces us to require that u o | • | lies 

in since we will apply Lemma El to the minimiser cjp on the right-hand side of E2), and we 

do not know any higher integrability property of this function than that (jP G 

It is clear that (imni remains true if W is replaced by kF — /, where / G Cb(M'^), since adding a 
constant to W — f does not change the value of the expression in the supremum on the right-hand 
side of EH), and the potential W — f — inf / also satisfies Assumption (W). 

Now we proceed with the proof of the lower bound of El) in Proposition El Go back to (imi) 

and recall that the random potential qj defined in (jd.d.SB is nonpositive and lies in Using 

Lemma El with h = qn and W replaced by VF — /, the last term on the right hand side of EH is 

estimated as follows. For any cf) G n satisfying ||<^||2 = 1, 


]g(^) ^ ^_/3j^(02)+/3(g^_w+/,02^^ 

Using this in EH; we obtain 


E 




> E'^^ 


I ie(^-2 ) 


:^-HN-2,f3-Kisr-2,f3+P{f,I^N-2,/3) 




(3.50) 


(3.51) 


Now we apply the same reasoning to the last expectation and iterate the argument. In this way we 
derive 


1 


Since u o I • I lies in 


-logE 

NP ^ 


> -Jp{p^) - (IT - /, * P\ P^) 


(3.52) 


lim {N — l)FAr * p = p‘ 

N^OO 


by assumption and since p'^ G L^(M'^), we have that 

2 _ j ,2 / y{\^x\) dx = p^‘iT:a{v) strongly in L^(M'^) 


Hence, 


lim *p‘^,p^) =d-KOi{v)\\p\\\. 

N—^oo Z 


(3.53) 
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This implies that 


lim inf-loff E 

Nf3 ^ 




^ -Xa{v)if)^ 


and the proof of the lower bound in (EU) is finished. 


A7Ta{v)\\4>\\l 


(3.54) 

□ 


4. Appendix 


4.1 Large deviations. 

For the convenience of our reader, we repeat the notion of a large-deviation principle and of the most 
important facts that are used in the present paper. See |DZ98] for a comprehensive treatment of this 
theory. 

Let A denote a topological vector space. A lower semi-continuous function I: A —> [0, oo] is called a 
rate function if I is not identical oo and has compact level sets, i.e., if /“^([O, c]) = {x G A: I{x) < c} 
is compact for any c > 0. A sequence (A 7 v)AfGN of A-valued random variables satisfies the 
large-deviation upper hound with speed a^ and rate function I if, for any closed subset F of A, 

limsup — logP(AAr G F) < — inf /(x), (4.1) 

N^oo O.N a:GF 

and it satisfies the large-deviation lower bound if, for any open subset G of A, 

liminf — logE(Av G G) < - inf I(x). (4.2) 

N—>oo aj\i x£G 

If both, upper and lower bound, are satisfied, one says that (A 7 v)Ar satisfies a large-deviation principle. 
The principle is called weak if the upper bound in (EH) holds only for compact sets F. A weak principle 
can be strengthened to a full one by showing that the sequence of distributions of Av is exponentially 
tight, i.e., if for any M > 0 there is a compact subset Km of A such that P(Ajv G M^) < for 

any n G N. 

One of the most important conclusions from a large deviation principle is Varadhan’s Lemma, which 
says that, for any bounded and continuous function F: A —> M, 

lim ^rlog [ dF = — inf (/(x) —F(x)). 

All the above is usually stated for probability measures P only, but the notion easily extends to suh- 
probability measures P = Pat depending on N. Indeed, first observe that the situation is not changed 
if P depends on N, since a large deviation principle depends only on distributions. Furthermore, the 
connection between probability distributions Pv and sub-probability measures Pat is provided by the 
transformed measure PAr(AAr G A) = PAr(AAr G A)/Fj\f{Xf^ G A): If the measures Pat o Xfj^ satisfy a 
large deviation principle with rate function I, then the probability measures Pat o Xf^^ satisfy a large 
deviation principle with rate function I — inf I. 

One standard situation in which a large deviation principle holds is the case where P is a proba¬ 
bility measure, and Xn = -^{Yi Yv) is the mean of N i.i.d. A-valued random variables Yi 

whose moment generating function M (F) = f e^Ai) jp ig finite for all elements F of the topological 
dual space A* of A. In this case, the abstract Cramer Theorem provides a weak large deviation 
principle for (AAr)AfGN with rate function equal to the Legendre-Fenchel transform of logM, i.e., 
I{x) = SUPA’e;t’*(F(x) - logM(F)). 
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4.2 Gross-Pitaevskii theory. 

Consider the ground-state energy per particle of the Hamilton operator H]\f in dm, 

1 


Xn = — inf 


2 = 1 


{h,HNh). 


(4.3) 


It is standard to show the existence, uniqueness and some regularity properties of the minimiser 
/iTv £ The large-behaviour, in a certain dilute regime, of xn and of the minimiser Hn was 

studied by Lieb et al. in a series of papers pYool . pm] . PyoT| . pp . see also the monograph 
[LSSY05] . It turned out there that the Gross-Pitaevskii formula in well approximates the 

ground-state energy. A summary of the large-A" results for xn is as follows. Assume that d G {2,3}, 
that u > 0 with u(0) G (0, oo], and v{r)r^~^ dr < oo, where a = infjr > 0: v{r) < ooj G [0, oo). 
These assumptions guarantee that the scattering length, denoted by 5(u), is finite l ILSSYOBj l. Note, 
that a{v) > a{v) ( jABKOd] 1. 


Theorem 4.1 (Large-A asymptotic of xn in 4 G {2,3}, ILSYODj . |LYnij . |LSYnij l. Replace v by 
vn{') = ■/9)v^) with = 1/A in d = 3 and 13‘j^ = A||(/>^|^j || in d = 2. Define 

G as the normalised first marginal of h‘j,f, i.e., 

= / h %{ x , X2 ,... , XN ) dx2 --- dxN , X G E'^. 

JlRd(iV-l) 


Then 


lim 

N^OO 


(GP) 

XN = 


and (jfijq *11 weak L^(Efi)-sense. 


The proofs show that the ground state, hjv, approaches the product-state if A gets large. 

Moreover, on the basis of this result, the occurrence of Bose-Einstein condensation in the ground- 
state (zero-temperature) was proved in pp . 
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